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ABSTRACT 

We calculate the joint period-spindown (P — P) distributions of millisecond radio 
pulsars (MSRP) for the standard evolutionary model in order to test whether the ob- 
served MSRPs are the unequivocal descendants of millisecond X-ray pulsars (MSXP). 
The P — P densities implied by the standard evolutionary model compared with ob- 
servations suggest that there is a statistically significant overabundance of young/high 
magnetic field MSRPs. Taking biases due to observational selection effects into account, 
it is unlikely that MSRPs have evolved from a single coherent progenitor population 
that loses energy via magnetic dipole radiation after the onset of radio emission. By 
producing the P — P probability map, we show with more than 95% confidence that 
the fastest spinning millisecond pulsars with high magnetic fields, e.g. PSR B1937+21, 
cannot be produced by the observed MSXPs within the framework of the standard 
model. 

Subject headings: X-rays: binaries — stars: neutron — stars: statistics — pulsars: 
general — pulsars: individual (B1937+21) 



1. Introduction 

Millisecond pulsars are commonly believed to be descendants of normal neutron stars that 
have been spun-up and recycled back as radio pulsars by acq uiring angular mom entum from their 
companion during: th e low-mass X-rav binary (XM XB) phase (A koar et al.lll98 2; Radhakrishnan & 



Srinivasan 19821 ) 



There are about ~20 high confidence nuclear or accretion powered (see Table 1) millisecond X- 
ray pulsars (MSXPs) which are thought to be the progenitors of millisecond radio pulsars (MSRPs) 



(jWijnands &: van der Klislll998l ). These MSXPs may become observable in radio wavelengths once 
accretion ceases, or the column density of the plasma from the fossil disk around the neutron star 
becomes thin enough to allow vacuum gap formation that leads to the production of coherent 
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radio emission. Towards the end of the secular LMXB evolution, as accretion rates fall below a 
critical value above wh ich detection presumably may be hampered due to absorption or dispersion 
([Thompson et al.lll994l ). the neutron star can re-appear as a MSRP. 



Although the connection between LMXBs and MSRPs has been significantly strengthened 
after the discovery of quasi-periodic kHz oscillations and X-ray pulsations in some tran s ient X -ray 



sources (jWiinands &: van der Klisl Il998l ; iMarkwardt et al 



radio pulsations from MSXPs have been detected so far (jBurgav et al.l 12003 ) 
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where Bg and Rq are the neutron star surface magnetic dipole field and radius in units of 10 9 G 
and 10 6 cm respectively. The Eddington limited accretion rate M^dd f° r a neutron star typically is 
~ 10~ 8 Mq yr^ 1 above which the radiation pressure generated by accretion will stop the accretion 
flow. This equilibrium period combined with the dominant mechanism for energy loss delineates the 
subsequent kinematics of the spun- up millisecond pulsar. The magnetic dipole model then implies 



a "spin- up region" (P ~ Pq^) (see lArzoumanian et al. Ill 9991 } on which the recycled neutron stars 



will be reborn as MSRPs. At the end of the active phase, MSXPs accreting with m and spinning 
with P eq , presumably transition into a MSRPs with an initial spin period of Pq ~ P ( 



eq- 



In the standard spin-down model, the MSRP evolution is driven by pure magnetic dipole ra- 
diation, i.e. braking index n = 3 in vacuum (see Manchester fc Taylorl 119771 : lLorimer Kramer 



2004). Alternative energy loss mechanisms such as multipole radiation or gravitational wave emis- 
sion, especially dur ing the initial phases of the reborn millisecond pulsars, have been suggested 



by several authors (IKrolik 



1991 



Bildstenlll998l ) but have yet to be observationally corroborated. 



Advanced Laser Interferometer Gravitational Wave Observatory (LIGO) will be able to probe the 
frequency space at which millisecond pulsars are expected to radiate gravitational waves, thereby 
putting stringent constraints on the micro physics of millisecond pulsars. 

The advances in radio observations, increased sky coverage with deep exposures of current 
surveys combined with robust post-bayesian statistical techniques that incorporate minimal as- 
sumptions, give us unprecedented predictive power on the joint period-spindown (P — P) and 
implied magnetic field (B) distributions. 

In this Letter, we attempt to go beyond phenomenological arguments and test whether MSXPs 
can produce the characteristics of the o bserved MSRPs within the framework of the standard model 



(jBhattacharya &: van den Heuvellll99ll . and references therein). 
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2. The Joint Period - Spindown (P — P) Distribution 

The evolution of millisecond pulsars can be consistently described in terms of i) the equilibrium 
period distribution {D) of MSXPs at the end of the LMXB evolution ii) the mass accretion rates 
(M) of the progenitor population during the recycling process iii) Galactic birth rates (P), and 
iv) the dominant energy loss mechanism after the onset of radio emission. 



2.1. Statistics 

We devise a semi-analytical evolution function £ to parametrize the evolution of millisecond 
pulsars after the accretion phase, which can be described in closed form as: 



Y,£(D l ,M t ,R i \at^) ^VVF(P,P) (2) 

i=0 



where VDT is the probability distribution function. The shap e parameters a an d (3 define the 



distributions (i.e., D,M,R for k=l,2,3) for the Beta functions 1 ( Evans et al. 2000 ) inferred from 



The evolution function £ is built by randomly choosing initiation seeds from a period distribu- 
tion D, which is then convolved via the standard model to consequently sample the P— P parameter 
space. For the observed MSXPs, the period distribution which seeds will be randomly chosen from 
is the observed Pmsxp distribution (table 1). We uniquely construct a "relaxed multidimensional 
Kolmogo rov-Smirnov (K-S) filter" (fi g. [T]) to check population consistencies by calculating the 



2D K-S (jFasano Sz Franceschinil 119871 ) probabilities (P2DK-S) between observed MSRPs and the 
synthetic population that is formed by these properly evolved progenitor seeds. The filtering is 
reiterated for each realization to obtain synthetic populations with consistent distributions as: 

D{a j j(3 } ) I^ Di (3 ) 

M(al(3?)^M t (4) 

R(aim^Ri (5) 

which is then used to construct the VDT in Equation [2j 

Nominally any Pidk-s > 0.2 value would imply consistent populations in a 2D K-S test. 
By allowing 0.005 < P2DK-S < 0.2 with lower fractions (see fig. [T|), we oversample outliers to 
compensate for possible statistical fluctuations and contaminations. A peak sampling rate around 
the nominal acceptance value of P2DK—S ~ 0.2 is the most optimal scheme that prevents strong 



x Beta functions are commonly preferred in Bayesian statistics as the least restrictive and most flexible prior 
distributions. It can take the form of an uninformative (e.g. uniform) prior, a monotonic line, concave, convex, 
unimodal (e.g. normal) or any extreme combinations of these shapes. 
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biases due to over or under-sampling. The main goal for oversampling outliers and relaxing the K-S 
filter is to test whether the standard model can at least marginally produce very fast millisecond 
pulsars with relatively high magnetic fields like PSR B 1937+21. 

The predictive significance of the P — P distribution for the probability map (figEJ) is obtained 
from a Monte-Carlo run with r = 10 7 valid realizations that produce consistent synthetic samples. 
Whilst sampling the P — P space, no assumptions were made regarding the progenitor period 
distribution (D), the accretion (M), or the Galactic birth (R) rates. The filter (eq. O HJ © is 
implicitly driven by the observed MSRPs. 

Fig [2] shows the expected P — P distribution for the standard model assuming that MSRPs 
have evolved from a progenitor population similar to the observed MSXPs. We do not include 
MSRPs in globular clusters because the P — P values in these cases may not necessarily be the 
sole imprint of the binary evolution, but can be significantly changed by possible gravitational 
interactions due to the crowded field. To explore the extend of the effects of an unevenly sampled 
progenitor population, we also show the region in the P — P space that is sensitive to alternative 
Pmsxp distributions. The probability map is overlaid with the observed MSRPs. 

3. Discussion and Conclusions 

The discovery of millisecond pulsations from neutron stars in LMXBs has substantiated the 
theoretical prediction that links MSRPs and LMXBs. Since then, the recycling process that pro- 
duces MSRPs on a spin-up region from LMXBs, followed by spin-down due to dipole radiation has 
been conceived as the "standard evolution" of millisecond pulsars. However, the question whether 
all observed MSRPs could be produced within this framework has not been quantitatively addressed 
until now. 

The standard evolutionary process produces millisecond pulsars with periods (P) and spin- 
downs (P) that are not entirely independent. The possible P — P values that MSRPs can attain 
are jointly constrained by the equilibrium period distribution (D) of the progenitor population, the 
mass accretion rates (M) during the recycling process and the dominant energy loss mechanism 
after the onset of radio emission. 

In order to test whether the observed MSRPs can be reconciled with a single coherent pro- 
genitor population that evolves via magnetic dipole braking after the spin-up process, we have 
produced the predictive joint P — P distribution of MSRPs for the standard model. We did not put 
restrictions on any of the parameters that drive the evolution. Acceptable D, M and R values were 
implicitly filtered. We have relaxed the K-S filter (see fig. H]) in order to oversample outliers and see 
whether it is even remotely feasible to produce young millisecond pulsars, like PSR B1937+21 or 
J0218+4232, that have higher B fields. The color contours in Figure [2] represent the P — P densities 
for MSRPs that are direct descendants of observed MSXPs (i.e. initial spin periods Pq ~ Pmsxp)- 
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The standard evolutionary model is able to successfully produce the general demographics of 
older MSRPs. It fails, however, to predict the younger and fastest MSRP sub-population that have 
higher B fields. 

Accretion rates that MSRPs have experienced during their accretion phase deduced from 
observed P — P values, combined with the observed MSXP period distribution (D = Pmsxp) 
produces mostly older MSRPs, including MSRPs with spin-down ages r c > 10 10 yrs. Figure [2] 
shows clearly that the apparent enigma of millisecond pulsars with spin-down ages older than the 
age of the Galaxy is mainly a manifestation of very low accretion rates during the late stages of 
the LMXB evolution. 

On the other hand, no physically motivated Pmsxp distribution has been able to produce the 
whole MSRP population consistently. The observed period distribution of MSXPs is likely to be 
under-sampled due to observational selection effects. It is also possible that some neutron stars in 
LMXBs simply do not produce observable pulses. In order to understand how the predicted P — P 
distribution is affected by different MSXP period distributions, we have estimated the whole extend 
of the P — P region that is sensitive to the prior. The values that may be produced for different 
Pmsxp distributions are shown by the shaded areas in Figure [2j No MSXP period distribution 
could mimic the observed relative ratios of young/old pulsars with high B fields. The fraction of the 
observed young/old MSRPs with high B fields is higher than what the standard model predicts by 
several orders of magnitude. This may further be e xacerbated by strong selection effects that limit 



our ability to observe very fast millisecond pulsars (jHessels et al.l 120071 ). The choice of a standard 
K-S test instead of the relaxed 2D K-S only increases the statistical significance. Hence, we argue 
that young millisecond pulsars with higher magnetic fields (e.g. PSR B1937+21) are inconsistent 
with the standard model. 

Therefore, it is tempting to suggest that the fastest spinning millisecond pulsars, in particular 
PSR B1937+21, may originate from a different evolutionary channel. While it appears that ordinary 
magnetic-dipole spin down from a source population similar to the observed MSXPs is adequate to 
explain the great majority of observed MSRPs, the low final accretion rates that are required cannot 
be reconciled with the high accretion rates needed to produce the fastest, youngest pulsars. We 
believe that it is necessary to posit the existence of a separate class of progenitors, most likely with 
a different distribution of magnetic fields, accretion rates and equilibrium spin periods, presumably 
among the LMXBs that have not been revealed as MSXPs. Understanding this additional channel 
is clearly critical to develo ping a natural solution to the long lasting "birth rate problem" (see, e.g. 



Kulkarni & Naravanlll988l ). 



It is also possible that the standard evolutionary model fails at another point. For example, 
if MSRPs during some portion of their evolution lose energy through a dominant mechanism other 
than magnetic dipole radiation (e.g. multipole radiation, gravitational wave or neutrino emission), 
then the evolution of pulsars through the P — P diagram could be complex. 

A combination of the above mentioned factors (i.e. alternative progenitors and subsequent non- 
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standard radiation) are then likely to play a role in millisecond pulsar evolution. A MSXP period 
distribution that has sharp multimodal features coupled with non-standard energy loss mechanisms 
may be able to reconcile for the joint P — P distribution of millisecond pulsars. 

The re search presented here h as made use of the August 2008 version of the ATNF Pulsar 



Catalogue ([Manchester et al. \\200w . The authors acknowledge NSF grant AST-0506453. 
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Fig. 1. — The K-S probability distribution of the synthetic populations used to sample the joint P — 
P parameter space. The relaxed 2D K-S filter allows additional acceptance of populations with 0.05 < 
P2DK-S < 0.2 by oversampling the extreme outliers of the P — P distribution in order to probe for possible 
contaminations and extreme fluctuations. The distribution also shows how optimally the P — P parameter 
space is sampled with a peak sampling rate around the nominal acceptance value of Pidk-s ~ 0.2. The 
dotted line is the relaxed 2D K-S filter for the synthetic populations that is used to construct the predictive 
distribution in fig [2] The solid line is the conventional K-S filter that would only accept strictly consistent 
P — P samples. 



This preprint was prepared with the AAS IATgjX macros v5.2. 
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Table 1. Accretion and nuclear powered pulsars 



v spin [Hz] Pulsar 



619 4U 1608-52 

601 SAX J1750.8-2900 

598 IGR J00291+5934 

589 X 1743-29 

581 4U 1636-53 

567 X 1658-298 

549 Aql X-l 

524 KS 1731-260 

435 XTE J1751-305 

410 SAX J1748.9-2021 

401 SAX J1808.4-3658 

377 HETE J1900. 1-2455 

363 4U 1728-34 

330 4U 1702-429 

314 XTE J1814-338 

270 4U 1916-05 

191 XTE J1807.4-294 

185 XTE J0929-314 

45 EXO 0748-676 



Reference 



Hartma n et al. (2003) 
Kaaret et al. (2002) 
Markwardt et al. (2004 ) 
Strohmaver et al. (1997) 
Zhang et al. (1997) 
Wiinands et al. (2001) 
Zhang et al. (1998) 
Smith et al. (1997) 
Markwardt et al. (2002) 
Kaaret et al. (2003) 
Wiinands van der Klis (1998) 
Chakrabartv fc Morgan (1998) 
Kaaret et al. (2006) 
Strohm aver et al. (1996) 
Markwardt et al. (1999) 
Markwardt fc Swank (2003) 
Galloway et al. (20011 
Markwardt et al. (2003) 
Galloway et al. (20021 
Villarreal Strohmaver (2004) 



Note. — The millisecond pulsar progenitor seeds used to construct 
the cumulative synthetic MSRP population for the observed Pmsxp- 
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Fig. 2. — The P — P distribution of millisecond pulsars for the standard model vs. observed MSRPs. 
The color contours are the expected MSRP distribution for the observed MSXPs, with red representing 
the highest density. Any MSRP outside of the color contours cannot be produced by the observed MSXPs 
with more than 95% confidence. The dashed red line is the 68% confidence limit of the expected P — P 
distribution taking the observed MSXPs as the progenitor population. MSRPs within the shaded region 
may be produced by MSXPs with spin distributions different than what is observed. The area shaded with 
lines assumes a maximum spin frequency u max — 619 Hz, which the fastest spinning observed MSXP (i.e. 
4U 1608—52). The shaded area extends to the dotted region if the maximum s pin frequency for M S XPs i s 
allowed to be v max = 760 Hz which is the theoretical upper limit predicted by IChakrabartv et al.l (|2003l ). 
The blue and red filled circles are the observed MSRPs in single and binary systems respectively. The area 
outside of the shaded region is not sensitive to the prior, i.e. the observed MSRPs outside of the shaded 
area cannot be produced consistently by the standard model for any MSXP spin distribution with more 
than 95 % confidence. The spin-down values for the observed MSRPs are corrected for secular acceleration 



dShklovskiilll970l : ICamilo et al. |[l994h . The spin-up line (P ~ rhP^ 3 ) for ra = MEdd and the characteristic 



age line for r c = 10 yrs are shown with dashed and dash-dotted lines. 



